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Abstract. Motivated by the possibility that different versions of the laws of physics conld 
be realized within other nniverses, this paper delineates the galactic strnctnre parameters 
that allow for habitable planets and revisits constraints on the amplitnde Q of the primordial 
density flnctnations. Previons work indicates that large valnes of Q lead to galaxies so dense 
that planetary orbits cannot snrvive long enongh for life to develop. Small valnes of Q lead 
to delayed star formation, loosely bonnd galaxies, and compromised heavy element retention. 
This work generalizes previons treatments in the following directions: [A] We consider models 
for the internal strnctnre of the galaxies, inclnding a range of stellar densities, and find the 
fraction of the resnlting galactic real estate that allows for stable, long-lived planetary orbits. 
[B] For high velocity enconnters, we perform a large ensemble of nnmerical simnlations to 
estimate cross sections for the disrnption of planetary orbits dne to interactions with passing 
stars. [C] We consider the backgronnd radiation fields prodnced by the galaxies: If a galaxy is 
too compact, the night sky seen from a potentially habitable planet can provide more power 
than the host star. [D] One conseqnence of intense galactic backgronnd radiation fields is that 
some portion of the galaxy, denoted as the Galactic Habitable Zone, will provide the right flnx 
levels to snpport habitable planets for essentially any planetary orbit inclnding freely floating 
bodies (bnt exclnding close-in planets). As the valne of Q increases, the fraction of stars in a 
galaxy that allow for (traditional) habitable planets decreases dne to both orbital disrnption 
and the intense backgronnd radiation. However, the onter parts of the galaxy always allow 
for habitable planets, so that the valne of Q does not have a well-defined npper limit (dne 
to scattering or radiation constraints). Moreover, some Galactic Habitable Zones are large 
enongh to snpport more potentially habitable planets than the galaxies fonnd in onr nniverse. 
These resnlts snggest that the possibilities for habitability in other nniverses are somewhat 
more favorable and far more diverse than previonsly imagined. 
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1 Introduction 

One of the most important cosmological parameters is the amplitude Q of the primordial 
density fluctuations. In our universe, this quantity is small, Q ~ 10“^, as measured by 
the corresponding temperature fluctuations in the cosmic microwave background [1-3]. The 
parameter Q sets the initial conditions for the growth of structure in standard cosmological 
models, where dark matter provides most of the matter available for structure formation and 
dark energy provides most of the energy density of the universe [2]. 

A key question is why the amplitude Q has this particular small value. The temperature 
fluctuations observed in the cosmic background radiation are thought to arise from quantum 
fluctuations produced during the inflationary epoch [4] in the early universe. However, the 
observed smallness of the parameter Q requires a corresponding small parameter in the scalar 
field potential of inflationary models [5, 6]. Other choices of the scalar field potential generally 
tend to produce larger fluctuations and hence larger values of Q (see the discussion of [7]). 
Given the possible existence of other universes, perhaps as part of a multiverse, it is reasonable 
to expect that these other regions of space-time sample a distribution of values for Q. Previous 
authors [8, 9] have considered this possibility and placed bounds on the values of Q that allow 
for structure formation with favorable properties. In general, if the fluctuation amplitude Q 
is too small, then cosmological structure has difficulty forming, and the resulting objects 
could be too rarifled to cool and form stars [8, 10, 11]; in addition, they could be too loosely 


1 




bound to retain the heavy elements produced by stars. On the other hand, if Q is too large, 
then galaxies are much denser than those of our universe; they can be too dense to allow 
for planetary orbits to remain unperturbed over the long time scales necessary for biological 
evolution to run its course. 

The consideration of alternate universes — in this case with different values of Q and 
hence different types of galactic structures — is counterfactual in nature. The laws of physics 
are described by a collection of fundamental constants, but we currently have no dehnitive 
explanation for why these constants have their measured values. Similarly, the standard 
picture of cosmology requires a corresponding set of cosmological parameters, with specihc 
values required to match observations, but we have no a priori theory to predict these val¬ 
ues. One possible - but only partial - explanation is that our universe is part of a larger 
space-time structure, often called the multiverse [12-14], where each separate universe repre¬ 
sents a different realization of the laws of physics within the ensemble of possible universes. 
More specihcally, each separate region of space-time can have its own values for both the 
fundamental constants of physics and for the cosmological parameters. Our universe - our 
region of space-time - thus has one particular realization of these quantities, drawn from some 
underlying distribution of parameters. 

This possibility motivates the present paper: If the multiverse concept can provide even 
a partial explanation for why our universe has its observed properties, then we must hrst 
understand what range of parameter space allows for a universe in which life (more or less 
as we know it) can develop [9, 15-18]. Note that in order to determine the probability of life 
arising in a universe, we must not only determine the parameters that allow for life but also 
the probability of realizing those parameter values [19]; thus far, we have little understanding 
of the underlying probability distributions. As a result, this paper has the modest goal of 
constraining what range of galactic properties allows for the survival of habitable planets, 
and using the result to contrain the allowed values of the fluctuation amplitude Q. 

The question of structure in other universes, including constraints on the primordial 
fluctuation amplitudes [9], has been considered previously [8, 16, 20]. In general, these prior 
treatments adopt a global approach, where a large number of constraints are presented and 
calculations are carried out in an order-of-magnitude fashion (see also [21]). Building on these 
previous studies, this paper focuses on the specihc issue of how galactic environments can 
disrupt potentially habitable solar systems. This work explores the problem in greater detail 
by considering the internal structure of galaxies (including a range of densities), numerical 
simulations of the cross sections for orbital disruption, and the radiation helds provided by 
the background galaxy. This paper thus complements previous studies. 

Within our Galaxy, taking inventory of habitable planets has now become observationally 
possible. The hrst planet with an Earth-like mass that resides in the habitable zone of a main- 
sequence star has recently been detected [22]. Moreover, projections suggest that the fraction 
of all Sun-like stars that have Earth-like planets in habitable orbits is large, of order 10 
percent [23]. This fraction of potentially habitable planets could be even higher for smaller 
stars, which are far more numerous. Although long anticipated, a large population of Earth¬ 
like planets is rapidly becoming an observational reality. 

In universes with different choices of the primordial huctuation amplitude Q, the density 
of stars in the resulting galaxies varies. For large Q, galaxies collapse early, while the back¬ 
ground universe is dense; the resulting galaxies are more compact, so that the habitability 
of Earth-like planets can be compromised in two conceptually different ways: [A] Potentially 
habitable planets can be scattered out of their solar systems when the stellar density is suf- 
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ficiently high. [B] The planetary surfaces can become too hot due to the background stellar 
radiation fields. As an order of magnitude estimate, both of these mechanisms start to com¬ 
promise the survival of habitable planets when galaxies are denser than those of our universe 
by a factor of ~ 3 x 10^. However, galaxies have internal structure, so that the stellar density 
varies enormously within such an extended stellar system. As a result, we consider galactic 
models and find the fraction of the stars within the galaxies that allow for habitable solar 
systems. This generalization allows some fraction of solar systems to survive, even in galaxies 
with extreme properties. 

The rest of this paper is organized as follows. The current paradigm for cosmological 
structure formation is briefly outlined in Section 2, where we discuss how the properties of 
dark matter halos scale with the amplitude Q and then specify models for galactic struc¬ 
ture. Section 3 presents results for the disruption of solar systems by scattering encounters, 
including calculation of interaction cross sections, disruption rates, and the fraction of solar 
systems that survive. Section 4 presents the analogous results for disruption due to radiation 
by calculating the background radiation fields of the galaxies and determining the fraction 
of solar systems that avoid overheating. The background radiation fields of dense galaxies 
provide regions where freely floating planets have the proper temperature to support liquid 
water on their surfaces; the extent of these Galactic Habitable Zones is explored in Section 
5. Finally, we conclude in Section 6 with a summary of our results and a discussion of their 
implications. 


2 Galaxy Formation and Galactic Structure 


2.1 Formation of Galaxies 

In this section, we use the basic cosmological framework for structure formation, as developed 
for our universe, to determine the types of structure that could form in other universes with 
alternate values of the fluctuation amplitude Q. Since we want to consider cases that are 
markedly different from our own, we use only the most basic elements of the current theory. 
To start, we assume that the universe is spatially flat and that the matter content is dominated 
by collisionless dark matter. In this context, the basic elements of structure formation with 
varying Q have been outlined previously for universes containing only dark matter [8] and 
including dark energy [9]. The derivation presented below follows these earlier papers (see 
[8, 9] for further detail). 

Density fluctuations begin to grow at the end of the radiation dominated era and become 
nonlinear after the overdensity reaches a critical value. For the case of a top-hat density 
perturbation with spherical symmetry, linear perturbation theory [24] predicts that nonlinear 
collapse ensues at an overdensity of ~ 1.7 (see also [25]). At the epoch of matter/radiation 
equality, the temperature Teq of the universe is given by 


Teq = rjmp— , 

i iv- 


( 2 . 1 ) 


where mp is the proton mass, rj is the baryon-to-photon ratio, Qm is the dark matter energy 
density relative to the critical density, and Hb is the corresponding energy density in baryons 
[9, 26]. Note that we are working in units where c = h = k = 1, so that G = (i-e., Mpi 

is the Planck mass). Notice also that Teq is essentially the mass in non-relativistic matter per 
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photon J The age of the universe at equality is then given by 


l/sy-'" Mpi _ Mpiiig 

~ 8ay(,™pSi„)^ ■ 

The mass scale of the horizon at the epoch of matter domination plays an important role and 
is given by 

Meq ^ « 4 X lO^^GeV 4 x , (2.3) 

where we have used t] = 10“® and = 6, values appropriate for our universe, to obtain 

the numerical estimates. Fluctuations cannot grow before the epoch of matter domination, 
so that all mass scales M < Mgq start with roughly equivalent fluctuation amplitudes and 
grow to become virialized at the later time given by 


teqQ 


- 3/2 f . 
Jvir ) 


(2.4) 


where the dimensionless parameter /vir 1 is a slowly varying function of the halo mass [8]. 
In actuality, the virialization time (2.4) becomes independent of M only in the limit of small 
mass scales; by using this approximation, we are thus encapsulating the mass dependence 
of the virialization time into the parameter /vir. Galaxies with smaller masses collapse to 
somewhat higher densities because they collapse earlier, when the universe was denser, and 
the density of the background universe at the time of collapse is a controlling factor in the 
problem (see Section 2.2 and especially Figure 1 of reference [8] for further discussion). 

After a galactic halo structure collapses, it has a characteristic density given by a factor 
of fc ~ 187r^ times the density of the background universe at time Gir; he.. 

Pc = fcPiUir) = ISTT^eq = I87r\2a rT^^). (2.5) 

\ ^vir / 

Inserting typical values, we obtain 

Pc ~ 5.5 X 10-14 g cm-3 gVv7' = 8-2 x IO^Mq pc'^ . (2.6) 


For mass scales M < Mgq, the characteristic densities of the collapsed structures are roughly 
comparable and given by equations (2.5) and (2.6). The density increases slowly with de¬ 
creasing galactic mass due to the factor /vir. Since we are interested in constraints produced 
by the densest galactic structures, and since galaxies in our universe have masses M < Meq, 
this paper focuses on this “low-mass” regime. 

To obtain numerical values for equation (2.6), we have used the baryon-to-photon ratio 
r] = 10“® and the mass density ratio ^M/^h = 6, although these quantities are expected to 
vary from universe to universe. As a result, one should keep in mind that the characteristic 
density scales according to the relation pc oc For ease of presentation, the results of 

this paper are given as constraints on value of Q, but the results more generally apply to 
the composite parameter QTgq^. In other words, larger values of the fluctuation amplitude 
Q can be offset by invoking alternate values of (p, f^b) that lead to smaller values of Teq 
(keeping the composite parameter constant). On the other hand, inflationary scenarios tend 

^Finally, note that we ignore the contribution of neutrinos in defining Teq. 
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to produce a wide range of possible fluctuation amplitudes, often larger than those in our 
universe, so it makes sense to use Q as the primary variable. 

For completeness, we note that for mass scales M > collapse occurs later, so that 
these large structures become nonlinear (as density perturbations) at later times given roughly 

by 

t ^ teqQ“"/Vvir(M/Meq) . (2.7) 

Alternatively, the mass scales of these larger structures are a fraction of the horizon mass, 

M ~ Q^/^Mhor ~ -^P'^hoT ) (2-8) 

where rhor is the horizon scale at the time when the perturbations become nonlinear. The 
corresponding length scale r is given by 

r ~ • (2.9) 

These more massive structures have lower characteristic densities by a factor of ~ (Meq/M)^. 
As a result, the lower-mass regime considered above produces galaxies that are more disruptive 
to their constituent solar systems. 

The considerations outlined above apply to the initial formation of cosmological struc¬ 
tures. As the universe in question evolves, structures on all scales tend to interact and merge, 
thereby creating a complicated network of interactions. Merging continues until the dark 
energy component of the universe (if nonzero) dominates and effectively freezes out further 
structure formation [9, 27-29]. While operative, the merging process tends to produce ever- 
larger structures and acts to redistribute the mass — both the internal structure of individual 
halos and the mass distribution of the ensemble of halos. In general, merging structures re¬ 
sult in more massive structures with lower density. Since we are primarily interested in the 
densest structures, however, we focus on halos with masses M ^ Meq. 

2.2 Structure of Galaxies 

Numerical simulations of cosmological structure formation indicate that dark matter halos 
of galaxies and clusters approach a nearly universal form. The benchmark study of this 
convergence [30] showed that the density distribution of the halos assume what is now called 
the NFW profile, which can be written as 

P = ? = —, ( 2 . 10 ) 

where rg is the scale length of the system. This simple form for the density distribution 
cannot continue out to arbitrarily large radii (as the enclosed mass would diverge). However, 
a number of additional effects act to provide halo edges. At the present cosmological epoch, 
halos with this density distribution meet up with their neighboring halos, which provide an 
effective outer boundary. At later times, however, the accelerating cosmic expansion leads to 
the isolation of halos, which approach an asymptotic form [28, 29] with a somewhat steeper 
density distribution with the form of a Hernquist [31] profile, 

P = ? = — • ( 2 . 11 ) 

These same simulations show that a truncation radius develops in above density distribution, 
where the asymptotic value vt ~ 4.6r2oo (see Figure 3 of [28]; see also [32]). 
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In this work, we assume that dark matter halos have the form given by equation (2.11). 
For simplicity, we also assume that the stellar component of the galaxies have the same form, 
albeit with a more concentrated configuration. In other words, we also use the Hernquist 
profile to specify the distibution of stars for a model galaxy. We can then calculate the 
density as a function of radius. At sufficiently large distances, the density will decrease 
enough so that planets can survive the disruptive effects of both scattering and radiation. 

The (crude) theory of structure formation outlined above shows that for a universe with 
a given value of Q, all halos with moderately low mass M < Meq will virialize at roughly 
similar times and have densities pc given by equation (2.5) after collapse. Here we make the 
identification that 

Pdm — Pc 1 (2.12) 


where pdm = Po is the density scale appearing in the density distribution (2.11) for the dark 
matter. For a given total halo mass M, the corresponding length scale ram for the halo is 
given by 


= 



MV'^ A-V 

Meq/ 87rr|q \9aRJ 


(2.13) 


With the parameters of the dark matter halo specified, we need to determine the corre¬ 
sponding parameters for the baryonic component. Since gas can dissipate energy, one expects 
the baryons to collapse further and attain a more centrally concentrated configuration. Be¬ 
cause the ratio of the total mass in baryons to that in dark matter is determined by 
only one additional parameter needs to be specified. We take this parameter to be the ratio 
Rp of the density scales of the two components, i.e.. 


Rp = 


Ph 

Pdm 


(2.14) 


Keep in mind that pb (respectively, pdm) is the density scale appearing in the baryonic (dark 
matter) density profile, both of which are assumed to have the Hernquist form given by 
equation (2.11). The scale length rb for the baryonic component is then given by 


rh = rdm 


Pdm\^^^ ( i4b 

-I = rdm I ^ 


Ph J 


K^uRf. 


1/3 


(2.15) 


If the star formation rate is uniform across the galaxy, the number density re*(^) of stars also 
has the form 

no 


n* = 


where ^ = — . 

rb 


The total number of stars in the (spherical) galaxy is given by 

W = 2TTrlno , 

where the stellar density scale no is given by 

Esf Ph 


no = 


(M,) ’ 


(2.16) 


(2.17) 


(2.18) 


where Cgf is the star formation efficiency of the galaxy (the fraction of the baryonic mass 
within the galaxy that has been processed into stars) and (M*) is the mean stellar mass. 
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This treatment for the baryonic component makes several assumptions: First we assume 
that the gas can cool and successfully make stars. In order to cool promptly, on time scales 
comparable to the age of the universe at the epoch of formation, the gas must be sufficiently 
dense. In this context, previous work shows that for Q ^ 3 x 10“®, the gas can cool for 
essentially all mass scales of interest [8]. As a result, the cooling criterion is expected to be 
satisfied for most galaxies in the large-Q universes considered in this paper. 

In our universe, during galaxy formation the gas often falls inward to form disk-like 
structures, rather than the simple spherical profiles considered here. In addition to providing 
a useful working model for galactic structure, the quasi-spherical density profiles invoked 
here are likely to arise for a number of reasons: The characteristic length scales for galactic 
structure decrease with increasing Q. In our galaxy, and most other spirals in our universe, the 
inner regions (on small length scales) display nearly spherical geometry in the form of galactic 
bulges. In fact, the Hernquist profile (equation [2.11]), which characterizes the asymptotic 
form of dark matter halos [28], was originally developed as a model for galactic bulges [31]. 
This form for the density profile also arises in a wide variety of other cold collapse scenarios 
(e.g., [33, 34]). If cooling and star formation proceeds rapidly, galaxies are expected to 
resemble scaled-up bulges, rather than disks. Even if they form initially, disk structures 
can relax dynamically by scattering stars and adjusting their structure into more rounded 
configurations [35]. For our galaxy, the disk will relax over a long time scale of order ~ 10^® yr 
[36, 37]. For the densest galaxies under consideration, the relaxation time can become shorter 
than the typical time for habitability (usually taken to be ~ 1 Gyr), but will often be longer. 
However, a particular kind of orbit instability — namely resonant motion perpendicular to 
the galactic disk — can force stars (or parcels of gas) to leave their orbital planes in even 
mildly triaxial systems [38]. In many cases, the time scale for this type of instability to 
alter the phase space of a stellar system is much shorter than the time scale for two-body 
relaxation [39]. Finally, merging of smaller structures plays an important role in determining 
the form of bulges and ellipticals in our universe, and mergers could be even more important 
in the high-Q universes of interest here. This paper thus assumes that the combination of 
rapid collapse, mergers, orbit instability, and dynamical relaxation leads to rounded galactic 
structures. 

3 Scattering Constraints 

In this section we estimate the disruption probability for solar systems in the idealized galactic 
models described above. The rate T of scattering encounters between a solar system and 
passing field stars is given by 

F = n^av , (3.1) 

where n* is the stellar density, a is the interaction cross section, and v is the relative speed 
between the solar system in question and passing stars. 

We start with an order of magnitude estimate: In our Galaxy today [40], the stellar 
density n* ~ 0.2 pc“^ and the random encounter speed u ~ 40 km s“^ = 40 pc/Myr. At 
these speeds, the cross section for disrupting the orbit of Earth is approximately cr ~ 50 AU^ 
= 1.25 xl0“® pc^ [41], where disruption corresponds to an increase in eccentricity to values 
above e > 0.5 (with this eccentricity, the stellar flux varies by a factor of 9, nearly an order 
of magnitude, over the orbit). With these values, the benchmark interaction rate for the 
disruption of our Solar System becomes 

F = 10"®Myr-^ = 10"®Gyr-i . (3.2) 
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Although the time scale required for life to develop is not known, 4.6 Gyr is apparently 
long enough and the value of tc = 1 Gyr is often used as a standard time scale required for 
habitability [42, 43]. Adopting this value, we would have to increase the stellar density by a 
factor of ~ 10® in order for the interaction rate to be high enough to compromise habitability 
(keeping other parameters fixed). 

3.1 Interaction Cross Sections 

In this work, we need to calculate the cross sections for disruption of planetary orbits due 
to interactions between a solar system and passing stars. Previous work has explored these 
cross sections for the lower speeds appropriate for star forming cluster environments [44, 45] 
and for field stars in our galaxy [41]. For the dense, highly interactive galaxies of interest in 
this work, we need the disruption cross sections for higher encounter speeds. 

For the sake of definiteness, we focus on the case of a single planet system interacting 
with a single passing star^^ and take the system to be an analog of the Earth-Sun system: The 
host star mass Mhost = lAf©, the planet mass Mp = IM^, and the initial planetary orbit is 
circular (e = 0) with semimajor axis a = 1 AU. Even in this simplified setting, the interactions 
must be described by a large number of additional variables, including the mass M* of the 
passing star, the impact parameter b of the encounter, the relative velocity Vr of the two stars 
at infinite separation, two angles {9, (f>) that specify the orientation of the solar system with 
respect to velocity vector of the incoming star, and the phase angle of the planetary orbit at 
the start of the interaction. We approach this problem using a Monte Garlo scheme to sample 
the distributions of all of the aforementioned variables. The stellar masses of the passing stars 
are drawn from the initial mass function, which is assumed to have the same form as that 
found in our universe. The angles are drawn randomly such that the passing star is equally 
likely to approach from any direction and the phase angle of the orbit is uniform-random. 
The relative velocity is drawn from a Maxwell-Bolzmann type distribution of the form df /dv 
oc exp[—n^/n^]. The velocity scale Ub of the distribution characterizes velocity dispersion 
of the background galaxy, where the mean value {v) = 2nb/\/^ ~ 1.128nb. Note that the 
galaxy could, in principle, have additional rotational velocities. 

With the variables chosen for a given interaction, the Newtonian equations of motion 
are integrated numerically using a Bulirsch-Stoer algorithm. Note that the integrations must 
be carried out to high accuracy: We need to resolve the orbital elements of the planet at the 
end of the simulation; the planet carries energy of order whereas the energy of the 

interaction is much larger, of order The orbital speed Uorb ~ 30 km s“^, but we need 

to consider relative speeds as large as Ub ~ 10^ km s“^; the ratio of these energies is thus 
~ 10 “ 10 , gQ numerical scheme must conserve energy at an even greater level of precision. 
For these simulations, the time step is reduced accordingly so that this criterion is met (see 
also [41, 44, 45] for further discussion). 

For each given velocity dispersion of the background galaxy, we have performed an 
ensemble of numerical simulations of interactions between single planet solar systems and 
passing stars. The results are then used to find the cross sections for a given change in the 
orbital elements of the Earth-analog planet [41, 44, 45]. This paper reports the results from 
scattering interactions using five choices of the velocity Ub that sets the scale of the velocity 
distribution; specifically we use log^Q(nb/lkm s“^) = 2, 2.5, 3, 3.5, and 4 (i.e., the velocity 
scale Ub is distributed evenly in the logarithm from 100 to 10^ km s“^). The total number 

^For the dense galaxies of interest, binary systems act like two separate single stars because the encounter 
speeds are much larger than binary orbit speeds [45]. 
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Figure 1. Cross sections for the disruption of a single-planet solar system due to passing stars. 
Cross sections are shown for increasing the orbital eccentricity to values of 0.5 - 1.0, or greater, from 
top to bottom, where post-encounter values e > 1 correspond to planetary ejection. The dashed 
curves shows a simple analytic fit to the cross section for obtaining e > 0.5 (see text). 


of numerical experiments was of order one million, with a greater number of realizations 
required for larger velocities. With the simulations completed, we then post-process the 
results to determine the probabilities and the corresponding cross sections for any particular 
outcome to occur (details of this procedure are given in [41, 44, 45]). 

The resulting cross sections are shown in Figure 1 for galactic velocity dispersions in 
the range -Ub = 10^ — 10“^ km s“^. Even with the large number of realizations, the Monte 
Carlo sampling errors for these cross sections are about 5%. We also note that these cross 
sections are in good agreement with extrapolations from previous simulations performed at 
lower impact speeds [41, 45]. The figure shows the cross sections for increasing the eccentricity 
of the planet from its starting value (e = 0) to a post-encounter value of e = 0.5 - 1, where 
the e = 1 final states represent ejections of the planet (including both stripping the planet 
from its bound orbit and driving the eccentricity to unity so that the planet is accreted by its 
host star). For purposes of this paper, we consider e > 0.5 to represent sufficient disruption 
for compromising planetary habitability and thus adopt the upper curve in Figure 1 as the 
destruction cross section. 

Figure 1 shows that the disruption cross sections are steeply decreasing functions of the 
velocity dispersion (set by Vb) of the galaxy. The computed cross section for increasing the 
eccentricity to e > 0.5 (the case of interest) is comparable to the geometric cross section of 
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the orbit (^4 = ~ 3.1 AU^) for velocity dispersion n ~ wq = 1000 km s“^. As a result, to 

a good approximation, the cross section for disruption (Tdis can be fit with a simple function 
of the form 


O'dis 


o'o 

u{l + u) 


where 


V V 

Vo 1000 km s“^ ’ 


(3.3) 


where (Tq = 1 AU^. This expression captures the numerically discovered behavior that the 
cross section falls more steeply when the velocity dispersion v ^ vq = 1000 km s“^ and 
approaches the form a oc v~^ in the limit n —?■ oo. This latter asymptotic form can be derived 
analytically using the impulse approximation and assuming that the passing star acts on the 
planet alone (so that the interaction with the host star is negligible). This latter assumption 
only holds in the limit where the cross section is much smaller than the geometric cross section 
and hence in the high speed limit. 


3.2 Disruption Rates 

With the interaction cross sections and the structure of the galaxy specified, we consider 
the disruption rates for solar systems residing at varying locations. Within the context of 
this galactic model, some portion of the the inner galaxy will always be too dense for solar 
systems to survive (because n* oc 1/^). Similarly, some solar systems can survive in any galaxy 
provided that the system resides at a sufficiently distant orbit. The goal is thus to determine 
that fraction of the stellar population that could in princple harbor habitable planets. 

The disruption rate for solar systems is given as a function of position within the galaxy 
by the expression 

^ noav 

and the corresponding requirement that solar systems survive for a given time tc takes the 
form 

no<T?;tc < ?(l + 0^> (3-5) 

where ^ = r/rb. As outlined above, this paper adopts the characteristic time scale tc = 1 
Gyr [42, 43], which is a resonable estimate for the time required for biological evolution to 
occur. The reader can scale the results for different choices of this time scale. Using the 
fitting function for the interaction cross section from equation [3.3], this constraint takes the 
form 

raoo-oWA < ■^(1 + 0^(1 + ^‘)) (3-6) 

where the scaled speed tt is a function of galactic position. To evaluate the speed u, we use 
circular orbits (of the solar system, around the galaxy) in the potential of a dark matter halo 
with the density distribution of equation (2.11), and find 



^ GM y/^ (grb/rd^)V^ 

Vr-dmi^oV l + erb/rdn,- 

The factors of rb/rdm arise because the orbital speed is determined by all of the mass in the 
galaxy, which is dominated by dark matter, whereas the number density of stars (and hence 
^) depends on the baryonic component, which can have a different scale length. Here M 
is the total mass of the galaxy. The interaction rate is given by the relative speeds of the 
solar system with respect to passing stars. The stars are expected to have random velocities 
superimposed on their overall orbital motion, where these random speeds are a fraction of 
the circular speed. Since both stars are moving, however, the relative speed is larger than the 
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random speed. For purposes of this paper, we assume these effects cancel and use equation 
(3.7) to determine the relative speeds. In the future, a more detailed treatment of the relative 
velocities should be developed, along with considerations of the orbits through the galaxy 
[46]. 


3.3 Survival Fractions for Scattering 

Equations (3.6) and (3.7) define the constraint that must be met in order for a solar system to 
survive over the time tc (taken here to be 1 Gyr). To evaluate this constraint for a given galaxy 
mass M, we need to determine the stellar density scale no, and length scales rdm and for 
the dark matter halo and the baryonic component, respectively. Equation (2.18) specifies no, 
where we take the star formation efficiency Cgf = 0.5, the mean stellar mass (M*) = 0.5Mq, 
and the ratio Rp = 10. This latter choice reflects the fact that baryons dissipate energy and 
can reach higher densities than the dark matter (and is roughly consistent with the relative 
densities in our galaxy [40]). We can then evaluate no by using equation (2.5) to specify pc, 
which depends on the fluctuation amplitude Q. The scale length rdm for the dark matter halo 
is then given by equation (2.13) and the corresponding scale for the baryonic component is 
given by equation (2.15). The remaining parameter that must be specified is the additional 
factor /vir that enhances the density (decreases the scale length). As reviewed in [8], /vir ~ 1 
for large halo masses near M = Mgq and decreases to /vir ~ 0.03 for a galactic mass halo 
M ~ 2 X IO^^Mq (where this latter value is close to the mass of the Local Group [47]). Here we 
model this additional collapse factor with a simple scaling law of the form /vir = (M/Mgq)^/^. 
Note that more sophisticated models can be used, where the density scale is written in terms 
of the halo concentration, which is then fit to numerical simulations [48]. 

In a galaxy with a given total mass M, a fraction of the constituent solar systems will 
reside at sufficiently large galactocentric radii to survive disruption of their planet orbits. 
Figure 2 shows the resulting survival fractions for galaxies with masses in the range M = 
10^^ — IO^^Mq (from lower left to top right in the figure). The survival fraction is shown as 
a function of the fluctuation amplitude Q. As the value of Q increases, large scale structure 
forms earlier and the resulting structures are denser (characteristic of the background density 
of the universe at their formation epoch). Although solar systems are more easily disrupted 
for galaxies in universes with large Q, some fraction of the systems always survive. Note that 
we present results for density fluctuation amplitudes up to Q —)• 1, but the formalism breaks 
down before reaching this limiting value (see the discussion of Section 6.2; see also [8]). For 
completeness, we also note that the number of surviving solar systems in a galaxy of mass M 
is given by 

Asurvive — RsurVne e,iM/{M,), (3.8) 

where Cgf is the star formation efficiency and (M*) is the mean stellar mass. 

The survival fraction depends on galactic mass. Within the framework used here, galactic 
structures with smaller masses (e.g., M = 10^*^ Mq, lower left curve) collapse to a greater 
degree, have smaller values of /vir, and become more dense (see equation [2.5]). Even in the 
most unfavorable case, with the smallest galactic mass and the limiting value Q —?• 1, the 
fraction of surviving solar systems is ~ 0.003, so that the galaxy is expected to have ~ 3 x 10® 
viable systems. For a galaxy with a mass comparable to that of the Milky Way, M = IO^^Mq 
( center curve), the fluctuation amplitude must be larger than about Q ^ 0.01 in order for 
more than half of the solar systems to be disrupted; in the limit Q —)• 1, about 3 percent of 
the systems survive, corresponding to ~ 3 x 10® viable solar systems per galaxy. 
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Figure 2. Fraction of solar systems that survive disruption from scattering encounters in galaxies 
of varying masses as a function of the fluctuation amplitude Q. In universes with larger Q, structure 
forms earlier, galactic structures are denser, and solar systems are more easily disrupted. The five 
curves shown correspond to galaxy masses M = 10^° — (from lower left to top right). 


4 Radiation Constraints 


In this section we find the radiation fields prodnced by the dense galaxies nnder consideration 
and estimate the fraction of solar systems that are disrnpted throngh overheating. For this 
assessment, we assnme that the interstellar medinm of the galaxy is optically thin to its 
stellar radiation. This condition is expected to hold becanse these dense galaxies are likely 
to experience efficient star formation. If enongh gas and dnst remain in the galaxy, however, 
the interstellar medinm conld be optically thick, and the material wonld become even hotter 
in order to transport ont the Inminosity of the galaxy. This treatment thns represents a lower 
limit to the disrnptive inflnence of the backgronnd radiation. On the other hand, we are 
implicitly assnming that the galaxy remains yonng enongh that the stars have not yet bnrned 
ont [36, 37]. 

To start, we consider a rongh estimate for the effects of radiation. The total flnx Fq 
dne to the backgronnd galaxy, at a particnlar location, is given by the snm over all of the 
constitntent stars. 


N 


Fg = Y1 


L, 


dvrr? ’ 

j=l 3 


(4.1) 
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where rj is the distance from the given location to the jth star and Lj is its InminosityJ 
To obtain an order of magnitnde estimate, we hrst assnme that the stars are distribnted 
spherically. The total number of stars N is then given by 


N = 



(4.2) 


where R is the size of the stellar system. We can replace the sum in equation (4.1) by an 
integral, i.e.. 


Fg 


rR 


47rr^(irn* 


(L) 

4yrr2 


n^{L)R , 


(4.3) 


where {L) is the appropriate average value of the stellar luminosity (see the discussion below). 
In order for the background flux of the Galaxy to exceed that of the host star, we need 


n^{L)R > 


L 


or 


9 {L) 

dvrro > 1, 

1j 


(4.4) 


where w is the orbital radius of the planet (assumed here to be in a circular orbit). We can 
eliminate R in favor of the number N of stars in the system. 


A 2 2/3 



1/3 


1j 


(4.5) 


In the Galaxy today, this quantity is about 3 x 10“® (for the choice {L)jL — 10; see below). 
The radiation flux from the background galaxy thus becomes too intense for habitability when 
the stellar density is increased by a factor of ~ 3 x 10^ (note that this density enhancement is 
roughly comparable to that needed for scattering encounters to compromise planetary orbits, 
as found in Section 3). 


4.1 Radiation Fields within Galaxies 

In this section we construct the expected radiation helds for galaxies with extended structure, 
where we use the results from Section 2.2 to specify the density prohles. For a solar system 
located at the center of the galaxy, the radiation held from the background has the form 


Fg = 


^TrrlfdC 


no 


(L) 


?(1 + 0^ 47rr§^2 


r 

noro{L) / 

Jo 




^(1 + 0^’ 


(4.6) 


where the scale length tq = rb is that of the baryonic component. Note that the integral is 
logarithmically divergent. As a result, for a solar system exactly at the galactic center, one 
must include a cutoff radius. 

In practice, however, solar systems will not reside at the center. For a solar system with 
dimensionless radial position ^ = a, the expression for the background hux takes the form 


Fg 


i^noML) 



idi 

(1 + 0^ 



dfi 

— 2^a/i ’ 


(4.7) 


^Here and throughout the paper, Fq is the magnitude of the total flux passing through the top of the 
atmosphere, and is evaluated at the location of the planet. Note that the flux of radiation outward through 
the galaxy is much smaller. 
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Figure 3. Fraction of solar systems that survive disruption from radiation in galaxies of varying 
masses as a function of the primordial fluctuation amplitude Q. As the value of Q increases, structure 
forms earlier, so that galactic structures are denser and galactic background radiation is more intense. 
Results are shown for {L)/L = 10 (solid curves) and 1 (dashed curves). For each case, the five curves 
correspond to (total) galaxy masses M = 10^° — Mq (from lower left to top right). 


where [x = cos 9 (and 9 is the polar angle in spherical coordinates). The angular integral can 
be evaluated to obtain 

1 

^ Jq (1 [f + + 2Ca] - log [f + - 2^a] } . (4.8) 

The integral can be evaluated so the flux has the form 


Fg 


noro{L) f _ 2 log a ' 
2(a2 - 1) [ a2 - 1 ■ 


(4.9) 


4.2 Survival Fractions for Radiation 

Next we use the results outlined above to estimate the fraction of solar systems that will 
reside in regions of the galaxies where the background radiation field does not compromise 
habitability. The requirement that the radiation from the host star, at the location of the 
planet, exceeds that of the background galaxy takes the form 

(a^ - 1) 1 - ^2°^^ >noro2TTZu^^, (4.10) 

— IJ L 
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where a is the dimensionless radial location within the galaxy. 

To evaluate this constraint, we need to specify the quantities on the right hand side 
of equation (4.10). For this estimate, the number density scale no and length scale tq were 
determined as described in Section 3, where we considered survival due to planetary scattering. 
We also need to specify the relative brightness of the stellar population through the ratio 
(L)/L, which measures the power output of the average star. In our present-day universe, the 
mass to light ratios for typical galaxies are of order 1 — 10 in solar units [35, 40], where this 
accounting includes all contributions to the mass (i.e., dark matter, gas, and stellar remnants 
in addition to main-sequence stars). We thus expect the ratio {L) jL to be of order unity, 
or perhaps somewhat larger. Here we present results for (L)/L = 1 and 10. Note that the 
corresponding light to mass ratio for zero-age main-sequence stars would imply a much larger 
ratio {L)!L ~ 1000 [49]; the galactic average value, appropriate in this context, is lower 
because the high mass stars die out quickly and stellar remnants are included. We also note 
that the era of nuclear burning stars lasts for trillions of years [37], but not forever, so that 
the ratio {L)/L will eventually become much smaller. 

Figure 3 shows the resulting fraction of solar systems that survive the intense radiation 
helds of their host galaxies, shown here as a function of the amplitude Q of the primordial 
density fluctuations. These results are qualitatively similar to those obtained previously for 
disruption due to scattering of planetary orbits (compare with Figure 2). However, scattering 
interactions are somewhat more disruptive, so that the survival fractions shown in Figure 3 
(considering the radiation helds) are higher. The most extreme case shown corresponds to 
the smallest galaxy with mass M = 10^^ TTq, the larger value of the luminosity ratio (L)/L, 
and the limiting value of the huctuation amplitude Q —>■ 1; even in this case, however, 1% of 
the solar systems are projected to survive the disruptive effects of radiation. 


4.3 Allowed Regions of Galactic Parameter Space 

The previous sections outline the requirements necessary for a solar system to survive disrup¬ 
tion due to scattering by background stars (Section 3.3) and intense radiation helds (Section 
4.2). These results show the fraction of surviving solar systems as a function of the huctuation 
amplitude Q (see Figures 2 and 3). In this section, we delineate the regions within a galaxy 
that allow solar systems to survive these two channels of disruption. 

Survival from scattering requires that the scattering optical depth ts is less than unity. 


TS 


riQavtc 

w+w 


(4.11) 


Survival from radiation requires that the radiation received from the background galaxy is 
less than that of the host star, i.e.. 


_ ^ _ norpl-Kw'^ (L) f _ 21og^ 
~ ~ C'^-1 L [ 


(4.12) 


Because planetary orbits can have different semimajor axes, this condition must be considered 
as approximate. Moreover, as considered in Section 5, the background galaxy could provide 
all of the power for a habitable planet. The constraint of equation (4.12) thus applies to 
conventional habitability. 

Both the scattering optical depth ts and the radiation ratio TZ are functions of radial 
position ^ within the galaxy. These functions are shown in Figure 4 for varying choices of 
the density scale uq. The lower curves in Figure 4 correspond to values comparable to those 
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Figure 4. Disruption levels as a function of position within an idealized galaxy. The solid curves 
show the probablity of scattering an Earth-like planet (over tc = 1 Gyr) for galaxy models with varying 
density scales ng. The dashed curves show the ratio of the radiation flux from the background galaxy 
to that provided by a central sun-like star. For both cases, curves are shown for galaxies with density 
scales no = 1 (bottom), 10^, 10"^, 10®, and 10® (top), in units of pc“®. 


of the solar neighborhood in the galaxy today. The remaining curves correspond to denser 
galactic environments, where the stellar density parameter is increased by factors up to 10®. 

A number of trends are clear from Figure 4. First, for density enhancements of order 
10®, the range of allowed parameter space starts to be compromised. Second, disruption by 
scattering interactions is usually more important than disruption by radiation. However, in 
the outer parts of the galaxy, the radiation fields provide the dominant channel of disruption. 
This behavior arises because the radiation flux decreases with radius as 1/^^ whereas scatter¬ 
ing disruption varies as 1/^^ (where both results are evaluated in the limit ^ 1:^ 1). Finally, 
we note that the outer parts of the galaxy are always diffuse enough to allow survival in the 
face of both scattering and radiation, i.e., some fraction of the solar systems always survive. 

Within the context of the formalism used here, the outskirts of the galaxy allow survival. 
In practice, however, the galaxy does not continue outward forever, but rather must have a 
most distant star. The fraction of stars contained within the dimensionless radius ^ is given 

by 


NjQ _ e 1+2g 

N (1 + 0^ (1 + 0"’ 


(4.13) 


where N is the total number of stars. In order of magnitude, the outermost star in the galaxy 
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resides at a radial location ^ given by A^(^) = iV — 1, or eqnivalently, 


l + 2e _ 1 
(1 + 0^ “ iv ’ 


(4.14) 


which has the solntion 


^ = N-1 + y/{N - 1)2 + AT - 1 Ri 2A^. (4.15) 

In order for the galaxy to be so dense that even the ontermost star experiences disrnption 
throngh solar system scattering, the reqnirement becomes 

noavtc > (21V) [1 + 2Nf ^ {2N)^ . (4.16) 

For a typical galaxy, the right hand side of this expression is ~ lO^'^, so that the density scale 
no mnst be larger than that of onr galaxy by a factor of ~ 10'^® in order for no solar systems to 
snrvive. Using simple scaling argnments, the density hqocQ^, so that the largest enhancement 
factor is of order 10^^ (<C lO'^*^), so that galaxies are never dense enongh to render all of their 
solar systems nninhabitable. 

Note that the preceding calcnlation assnmes that the stellar component of the galaxy 
does not have a well-defined onter edge. In practice, however, the onter parts of the galaxy 
can become too rarified to form stars [50], where the intergalactic radiation field prevents the 
formation of molecnles. As a resnlt, the ontermost portions of the galaxy wonld have to be 
popnlated with stars via dynamical relaxation. 

Before leaving this section, we note that scattering and radiation are ronghly comparable 
as disrnptive inflnences on habitable solar systems for galaxies in the present-day nniverse. 
We can nnderstand this apparent coincidence as follows: The scattering optical depth is given 
by the prodnct noavtc (from eqnation [4.11]), whereas the analogons prodnct that determines 
radiative disrnption is ~ no^oTTro^ (from eqnation [4.12]). Both of these qnantities thns 
depend linearly on the stellar density scale uq. In addition, the cross section a for scattering 
is ronghly comparable to the area of the planetary orbit vrro^, which determines how mnch 
radiation is received from the central star. The coincidence is that the total distance traveled 
by a solar system dnring the benchmark time scale {vtc) is comparable to the distance scale 
ro of the galactic density profile (i.e., the orbit time ro/v ~ tc)- More compact galaxies — 
and especially their central regions — have shorter orbit times and are more disrnptive via 
scattering (compared to radiation). 

5 Galactic Habitable Zones 

The intense radiation fields provided by compact galaxies introdnce another possible channel 
for habitable planets. If the galactic backgronnd radiation has the proper intensity, then 
any potentially habitable planet can have the right temperatnre to snpport liqnid water 
on its snrface (one of the nsnal reqnirements for habitability [51]). More specifically, the 
snrface temperatnre of a planet will have a minimnm valne, independent of its orbit within 
its host solar system (inclnding nnbonnd orbits). As expected, and as shown in Section 4, 
the backgronnd radiation flnx within a galaxy decreases with galactocentric distance, so that 
a spherical shell within the galaxy will provide the proper radiations levels. We denote this 
region as the Galactic Habitable Zone (GHZ). 
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Figure 5. The Galactic Habitable Zone is the portion of a galaxy where the background radiation 
has the right intensity to support biospheres, independent of the orbital elements of the planet. The 
top panel shows the inner and outer boundaries of the GHZ as a function of the composite parameter 
X that encapsulates the relevant galactic properites (see text). The bottom panel shows the fraction 
/ of the solar systems (per galaxy) that reside in the GHZ. If each star produces on average 10 planets, 
then galaxies with / > 0.1 (above the dashed line) will support more potentially habitable planets 
than galaxies in our universe. 


5.1 Extent of the Galactic Habitable Zone 

We start by delineating the bonndaries of the GHZ. Section 4 shows that the backgronnd 
galaxy provides a radiative flnx Fq given by eqnation (4.9). By eqnating this flnx to that 
reqnired for habitability, we obtain a constraint of the form 


Fg 


C 


Lrr 


© 


47rti7^ 


(5.1) 


where w is the semimajor axis of the Earth-analog orbit and the factor C is a dimensionless 
constant of order nnity. The reqnirement of liqnid water on the snrface of a planet depends 
on the planetary atmosphere and other properties. Here we nse the reqnirement that the 
backgronnd flnx of the galaxy Fg is within a factor of two of that provided by the Snn to 
the Earth, i.e., we allow C to vary over the range 0.5 < C < 2. Next we define a composite 
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parameter X that incorporates the galactic properties, 
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(5.2) 


The Galactic Habitable Zone then lies between radial locations 
reqnirement 
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2 log I 
^ 2-1 


C 

X’ 




that statify the 


(5.3) 


where C = 1/2 for the onter bonndary and C = 2 for the inner bonndary of the GHZ. 

The resnlting bonndaries of the GHZ are shown in the npper panel of Fignre 5 as a 
fnnction of the composite parameter X. The corresponding fraction / = Nh/N^ of solar 
systems that reside in the GHZ are shown in the lower panel of the fignre. The horizontal 
dashed line delineates the portion of parameter space for which 10 percent of the solar systems 
reside in the GHZ, which corresponds to valnes of the parameter X in the range X ~ 0.4 — 125. 
Note that nearly all of the planets in the GHZ are potentially habitable, by definition. In 
contrast, within ordinary galaxies (like onr own), only a fraction of planets have orbits at 
the right distance from their host stars to snpport biospheres. Moreover, this fraction is only 
abont 10 percent: The observed orbital spacing in solar systems is snch that only abont one 
star per solar system has the right temperatnre, and solar systems can often snpport ~ 10 
planets. This resnlt raises an interesting possibility: Galaxies with composite parameter 
X = 1 — 100 conld snpport more habitable planets than those in onr nniverse. 

The condition of having a robnst GHZ reqnires the parameter X to lie in the range 1 

- 100, which reqnires that the colnmn density of stars in the galaxy no^'o ~ 10® — 10^® pc“2 
(where we take {L)/L = 1). For comparison, typical colnmn densities for galaxies in onr 
nniverse are of order uqVq ~ 100 (e.g., for a galaxy with N = 10^^ stars and scale length vq 

— 12.5 kpc). As another benchmark, the snrface brightness of onr Galaxy is ~ lOOOL© pc”^ 
[40], which converts to a colnmn density in stars of order 100 - 1000 pc“2. As a resnlt, within 
the context of this formalism, the most favorable galaxies for potentially habitable planets 
mnst have colnmn densities enhanced (relative to onr nniverse) by factors of 10® — 10®. 

We wonld like to find the range of flnctnation amplitudes Q that prodnces galaxies 
with large GHZs. Unfortnnately, there is not a definitive one-to-one mapping between the 
constraints on galactic properties and analogons constraints on the density flnctnation param¬ 
eter Q: Universes with a given Q produce galaxies over a range of mass scales, and different 
galactic masses collapse to a range of densities and sizes. Nonetheless, we can make a rongh 
estimate nsing typical valnes for the galactic masses. Since we expect the density to scale 
with flnctnation amplitnde according to uq oc Q® and the length parameter to scale according 
to ro oc 1/Q, the colnmn density scales as Q^. The preferred flnctnation amplitnde Qopt for 
habitable planets is thns larger than that of onr nniverse by a factor of 10® — 10^, so that the 
optimnm amplitnde is expected to fall in the range 


Qopt = 0.01-0.1. (5.4) 

The above scaling argnment indicates that a snbstantial GHZ reqnires relatively large 
valnes of the amplitnde Q. In snch nniverses, however, several complications shonld be kept in 
mind: Large flnctnation amplitndes can lead to rapid strnctnre formation, and conld increase 
the nonlinear mass scale beyond that realized in onr Universe. If the component of dark 
energy — which acts to accelerate the nniverse — is snfhciently small, then galaxies can 
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readily merge with each other, and the first generation of galaxies would not survive long 
enough to support life. The question of habitability would then depend on the nature of the 
resulting merged galactic structures. In addition, the central regions of these dense galaxies 
are likely to support the formation of massive black holes, which would alter the galactic 
luminosity. 

Some authors have suggested that the probability of an inflationary universe producing 
a given amplitude for the density fluctuations is an increasing function of Q [7]. This result is 
consistent with the finding that small values of Q require carefully chosen small parameters 
in the inflationary potential [5, 6 ]. As a result, universes with larger values of Q, including 
the range indicated in equation (5.4), could be more common than our own. 

5.2 Properties and Consequences 

For the larger stellar densities expected within the GHZ, planets have a good chance of being 
scattered out of their orbits. Over the range of galactic masses and the amplitudes Qopt 
given in equation (5.4), anywhere from 3% to 97% of the planets in Earth-like orbits will 
experience disruption (see Figure 2). Planets on wider orbits will be disrupted to an even 
greater extent. As a result, the majority of the viable planet population in this scenario will 
be freely floating, rather than bound in orbit about their host stars. Although scattering 
encounters are expected to liberate planets in the GHZ, stellar collisions will be rare for time 
scales tc- At the high densities of the galactic centers, however, stellar collisions can more 
readily take place. The centers of these galaxies are thus favorable sites for the formation of 
large black holes [52]. 

Freely floating planets, and entire solar systems, can eventually leave the GHZ through 
accumulated perturbations from passing stars. The time required for this process of dynamical 
relaxation is given by 

('c; c;\ 

= STTAOm^n, ’ 

where v is the orbit speed and A = ln[ 62 /bi], where bi and 62 are the minimum and maximum 
possible impact parameters [35]. Using equation (5.5) to evaluate the relaxation time for 
galaxies with X = 1 — 100 and for locations within the GHZ, we find 

*“■“ ~ (i^) ^ ^ '5) 

where N^, is the total number of stars in the galaxy. Since the dimensionless radius ~ 1 in 
the GHZ, the relaxation time is much longer than the characteristic time scale tc = 1 Gyr 
required for habitability. As a result, most potentially habitable planets will stay within the 
GHZ long enough for biological evolution to take place. ^ 

In a related but different context, the concept of a Galactic Habitable Zone has been 
put forth for our Galaxy [53]. In that case, however, the conditions required to compromise 
habitability are much less extreme than those of this paper. Specifically, the habitable zone 
from [53] is based on the chemical makeup of the galaxy as a function of position (and 
requires the metallicity to be within a factor of ~ 2 of the value for our Solar System). Since 
the metallicity Z decreases with galactocentric distance, the outer parts of the galaxy have 
low Z and are less suited for life. The definition of the GHZ used in this paper is similar in 

^Note that the central regions of the galaxies, where ^ ^ 1, have much shorter relaxation times, so that 
rearrangement of stellar orbits is expected there. 
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spirit to the constraints of [54], who considered the disruption of planetary systems in rich star 
clusters. In this latter study, sufficiently rich stellar environments produce intense radiation 
fields that heat planet-forming disks above the temperature required for ices to condense; 
with no ices, the process of giant planet formation can be compromised. 

This present work focuses on the survival of planetary systems after they are made; 
future studies should also consider how environmental effects can hinder their formation. In 
particular, the background temperature within the GHZ allows for water to be in liquid form 
at atmospheric pressures, which implies that the water could be gaseous at the lower pressures 
characteristic of planet-forming disks [55]. In this case, the dust grains - and later the rocks 
- that provide the building blocks for planets will have little water, and terrestrial planets 
could end up dry. Note that this dilemma is a scaled up version of the same issue facing 
potentially habitable planets in our universe: In the inner solar system, where planets are 
warm enough to support biospheres, the rocky planets must be built out of raw material with 
relatively little water. As a result, alternate channels for the delivery of water to Earth, and 
hence other terrestrial planets, have been considered. If enough of the nebula remains intact 
after a terrestrial planet is formed, its atmosphere can be both massive and rich in hydrogen. 
In this setting, the atmospheric hydrogen can react with various oxides (such as FeO) in the 
still-liquid rocky planet and produce a significant amount of water [56, 57]. In any case, the 
delivery of water to potentially habitable planets in the GHZ should be explored further. 

For completeness, we note that the radiation fields found in the GHZ are much less 
intense than those required to make the night sky as bright as a stellar surface, the condition 
considered in Gibers’ paradox [58, 59]. In order for any given line of sight to intercept a 
stellar surface, the optical depth through the galaxy would have to exceed unity, so that 
noroTri?^ ~ 1 (compare with the analogous expressions for disruption from equations [4.11] 
and [4.12]). However, stellar radii are much smaller than planetary orbits (ru ~ 1 AU) 
or scattering cross section length scales {£ ~ ^/o^ ~ 1 AU). As a result, the column density 
noro of stars through the galaxy would have to be much larger than the values considered 
here, by a factor of (ro/i?*)^ ~ 4 x 10^, for the night sky to be as bright as the surface of the 
Sun.^^ This distinction is important: The incoming radiation has the energy characteristic of 
stellar photospheres (with temperatures T* ~ 3000 — 10,000 K), whereas the planets have 
much lower temperatures Tp ~ 273 — 373 K (and emit lower energy, infrared radiation). 
This temperature difference allows for a heat engine to operate and thereby perform the work 
necessary to run a biosphere. 

In order for a biosphere to operate, the heat generated by the planet must be able 
to radiate away. In our present-day universe, the background temperature of the galaxy 
(Tgal r\j 3—10 K) is much colder than the planetary temperature Tp, so that exhaust heat poses 
no problem. Since the background temperature of the galaxy is hotter in the GHZ, one might 
worry that a heat engine could have trouble operating. As outlined above, however, the galaxy 
is optically thin to its internal radiation, including both stellar radiation and that re-radiated 
by planets. As a result, the ultimate heat sink for the galaxy is cold extragalactic space (where 
T <C Tp). Planets can thus absorb optical light and freely re-radiate the energy as infrared 
radiation. This process generates large amounts of entropy, which allows for a decrease in 
the entropy in parts of the biosphere itself. In addition, planetary atmospheres raise the 
temperature of the planets further, and this increase also helps to radiate the necessary 
exhaust heat. For completeness, we also note that a biosphere does not function as a simple 

Tn other words, in the GHZ, the night sky is (approximately) as bright as the daytime sky on Earth, but 
not as bright as a stellar photosphere. 
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carnot heat engine, so that its efficiency is not given by the classical expression (e = l—Tc/Th). 
Instead, one mnst consider the Gibbs free energy of the system, inclnding all of the terms 
associated with the different particle components and their chemical potentials (for farther 
discnssion, see [42, 43]). 

6 Conclusion 

Extending previons work, this paper has shown that when galactic stellar systems become 
snfficiently dense relative to those in onr present-day nniverse, significant disrnption of solar 
systems can take place. Snch disrnption can render Earth-like planets nninhabitable throngh 
two different mechanisms: Scattering interactions can pertnrb planetary orbits before life has 
time to evolve and radiation from the backgronnd galaxy can dominate that of the host star 
and thereby make all planets too hot to snpport biospheres. Nonetheless, some fraction of 
the potentially habitable systems will snrvive. A snmmary of onr resnlts is given below (in 
Section 6.1), followed by a discnssion of their implications and indications for fnrther work 
(Section 6.2). 

6.1 Summary of Results 

We have performed a large ensemble of ~ 10® nnmerical simnlations to determine the inter¬ 
action cross sections for the disrnption of planetary orbits via enconnters with passing stars 
(see Fignre 1). The disrnption cross section of Earth-analog orbits can be fit with the fnnc- 
tion given by eqnation (3.3). As expected, this cross section approaches the asymptotic form 
a ~ 1/n^ in the high speed limit, where the limiting regime reqnires v ^ 1000 km s“^. 

The density of stars within a galaxy is a fnnction of position; here we nse the Hernqnist 
profile of eqnation (2.11) as a working model. For this class of profiles, both the scattering 
rate T and the backgronnd radiation flnx Fq diverge in the inner limit ^ 0. On the other 

hand, both T —)• 0 and Fq —)• 0 in the limit of large radii ^ 3> 1. Inclnsion of galactic strnctnre 
thns changes the problem from previons treatments, which implicitly assnme that galaxies 
are characterized by a single density. In every galaxy, some fraction of the solar systems 
will be disrnpted by scattering enconnters, whereas some other fraction will be disrnpted by 
backgronnd radiation fields (and some fraction will snrvive). 

The characteristic densities for galactic strnctnre (see eqnation [2.5]) increase with the 
amplitnde Q of the primordial density flnctnations. As a resnlt, the fraction of solar systems 
within a galaxy that snrvive disrnption - and thereby remain potentially habitable - is a 
decreasing fnnction of Q. This trend is shown in Fignre 2 for disrnption of planetary orbits 
by scattering enconnters and in Fignre 3 for disrnption by the backgronnd radiation field of 
the galaxy. Note that the snrvival fractions also depend on galactic mass (as shown in the 
fignres). The onter parts of the galaxies remain habitable for essentially any choice of the 
galactic strnctnre parameters (see Fignre 4). 

Becanse of the intense galactic backgronnd radiation fields, some portion of the galaxy 
can provide the right flnx levels to snpport habitable planets. Any planets residing in this 
region - denoted here as the Galactic Habitable Zone - can have snrface temperatnres com¬ 
patible with maintaining liqnid water. Althongh planets residing too close to their host stars 
will still be too hot, planets in all other orbits, inclnding freely floating bodies, can in principle 
snpport biospheres. Moreover, the GHZ is large enongh so that some galaxies in other nni- 
verses can snpport more potentially habitable planets (powered by the backgronnd radiation 
of the galaxy) than the galaxies in onr present-day nniverse (Fignre 5). 
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6.2 Discussion 


In approximate terms, disruption of habitable planets through both scattering and radiation 
becomes important when the background density of stars is a million times larger than that 
of the Solar Neighborhood, which is not atypical within galaxies of our present-day universe. 
Since galactic densities increase with the approximate scaling p oc (equation [2.5]), we 
naively expect the disruption of habitable planets to become important when Q is about 100 
times larger than in our universe, i.e., for Q ^ 10“^. This expectation is borne out in Figures 
2 and 3, although universes can remain habitable for even larger values of Q if we only require 
that a fraction of the solar systems remain viable. 

The results of this work indicate that habitable planets can survive in galaxies with a 
wide range of properties, and hence in a wide range of universes. At first glance, this result 
may seem to contradict previous work that places relatively tight limits on the fluctuation 
amplitude Q. In this current generalized treatment, however, more planets can survive for 
several reasons: [A] Galactic structures naturally include a wide range of stellar densities, with 
lower densities in the outer regimes that allow for habitable orbits to survive. [B] With the 
consideration of internal galactic structure, we can determine the fraction of solar systems 
that avoid disruption (rather than using an all-or-nothing approach). Because planets are 
thought to be common, and galaxies contain many solar systems, a galaxy (universe) can 
remain habitable with only a small fraction of its solar systems surviving (note that Figures 
2 and 3 are plotted with logarithmic scales). [C] Even planets that are stripped from their 
host stars can remain habitable if they reside within the Galactic Habitable Zone, i.e., the 
region of the galaxy where the galactic background radiation has the proper intensity to heat 
planets. 

The results of this paper are presented for density fluctuation amplitudes as large as 
Q = 1. However, the reader should keep in mind that the formalism breaks down before 
reaching such high values. For sufficiently large fluctuation amplitudes, structure formation 
is thought to take place violently, thereby leading to rampant black hole formation and other 
departures from the standard scenario. This regime of structure formation has not been 
well-studied via numerical simulations, but such work should be carried out in the future. 
Even for the regime of moderately higher values of the amplitude Q, numerical simulations 
of structure formation should be carried out in order to verify the assumption that dark 
matter halos approach a nearly universal form (assumed here to be NFW/Hernquist). We 
also need a better specification of what happens to the baryons in these high-Q scenarios. In 
this work, we assume that galaxies become scaled-up versions of the galactic bulges in our 
universe (with density profiles of equation [2.11]). If the galactic system has enough angular 
momentum, however, baryons can settle into a rotationally supported disk structure (like the 
spiral galaxies of our universe). The survival fractions and background radiation fields should 
also be determined for galaxies with disk geometries. 

This treatment considers universes with different amplitudes Q for the density fluctua¬ 
tions, and hence denser galactic structures, but holds fixed the other properties of the universe. 
In general, however, both the fundamental constants of physics and cosmological variables 
(in addition to Q) could be different in other universes. Such additional variations would 
change the predicted fraction of habitable solar systems. As one example, if the gravitational 
constant G is larger, then stars are brighter [60-62], and habitable planets must have larger 
orbits; as a result, planetary scattering could be enhanced. On the other hand, this change 
could be offset with an increase in the fine structure constant a, which would act to suppress 
nuclear fusion and would require closer orbits to achieve chemical reactions [20, 63]. Future 
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work should thus take into account coupled variations in all of the relevant parameters. In 
addition, the ‘constants’ of nature could also vary with time within a particular universe 
[64, 65]. Considering all of the possibilities is an enormous undertaking; the hrst step is to 
consider more limited possible variations as carried out here. 

Finally, we note that most considerations of alternate universes use the properties of 
our own universe as the baseline for comparison, often with an implicit assumption that our 
universe provides the best possible environment for the development of life. The results of 
this paper call this assumption into question: With the right parameters, galaxies can support 
extensive Galactic Habitable Zones, where up to 20% of the stars, and hence up to 20% of all 
planets, reside within a radiation held comparable to the one the Earth receives from our Sun. 
Such galaxies - residing in universes with the right properties - could support more habitable 
planets than our own. Within our universe. Earth is often considered as the optimal planet 
for supporting a biosphere, but other planets could be even more favorable [66]. This issue, 
which could be called super habit ability, should also be explored further. 
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